arXiv:nucl-th/0512024v2 25 Jan 2006 


v4-dependence of hadronization in nuclei 

H.P. Bloki’^.Q and L. Lapikas^'lll 
^Department of Physics and Astronomy, Vrije Universiteit, 
de Boelelaan 1081, 1081 HV Amsterdam, The Netherlands 
^Nationaal Instituut voor Kernfysica en Hoge-Energie Fysica (NIKHEF), 

P.O. Box fl882, 1009 DB Amsterdam, The Netherlands 
(Dated: February 7, 2008) 

The ^-dependence of models for the attenuation of hadron production in semi-inclusive deep- 
inelastic scattering on a nucleus is investigated for realistic matter distributions. It is shown that 
the dependence for a pure partonic (absorption) mechanism is more complicated than a simple 
(^1/3) behavior, commonly found when using rectangular or Gaussian distributions, but that the 
A-dependence may still be indicative for the dominant mechanism of hadronization. 
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The study of hadronization, the process that leads 
from partons produced in some elementary interaction to 
the hadrons observed experimentally, is of importance, 
both in its own right as a study of a non-perturbative 
QCD process, and in the interpretation of data from ex¬ 
periments that use outgoing hadrons as a tag. The end 
products of the hadronization process in free space are 
known from e“ annihilation, but very little is known 
about the space-time development of the process. One 
way to investigate this is to study the semi-inclusive pro¬ 
duction of hadrons in deep-inelastic scattering of elec¬ 
trons from a nucleus, where the nucleus is used as a 
length(time)-scale probe (see Refs. Hi) 

Even if hadronization is not yet quantitatively un¬ 
derstood, it is known that the following processes play 
a role in leptoproduction of hadrons in a nucleus (see 
also Fig. ^). After a quark in a nucleon is hit by 
the virtual photon, it loses energy by scattering from 
other quarks and radiating gluons, thus creating quark- 
antiquark pairs. After some time^ and corresponding 
length If colorless (pre)hadrons^ are formed. The av¬ 
erage value T/ of the formation length has been esti¬ 
mated 0 based on the Lund model to be typically in the 
range of 1-10 fm if the virtual-photon energy is in the 
range 5-30 GeV. Hence, Lf is comparable to the size of 
a nucleus. 

If the hadronization is fast, i.e., the hadrons are pro¬ 
duced inside the nucleus, they can be absorbed, which 
will show up as an ‘attenuation’ of the hadron yield. (In 
experiments one measures the ratio Ra of the yield on a 
nucleus with mass number A and the one on a deuteron.) 
If, on the other hand, the hadronization is stretched out 
over distances large compared to the size of a nucleus. 
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^ For a more detailed discussion of the concept of formation time 
etc. see Ref. |3- 

^ For the present discussion it is not needed to discriminate be¬ 
tween hadrons and prehadrons, so in the remainder we will just 
talk about hadrons (but see, e.g., Refs. HH) 



FIG. 1: Illustration of hadronization in a nucleus and nsed 
coordinates. The parton is produced at point A, while the 
hadron is formed at B. The distance If has a distribution 
with as average the formation length Lf. 

the relevant interactions will be partonic, involving the 
emission of gluons and quark-gluon interactions, which 
also changes the production of hadrons. 

Since hadronization, as a non-perturbative QCD pro¬ 
cess, cannot be calculated from first princi]de^ various 
models have been developed (see, e.g.. Refs. HIMIaiHIE 
0) to describe hadron production and attenuation in a 
nucleus. Some models focus on the partonic part, while 
others include or emphasize the hadronic part. In all 
cases a sizable dependence on the mass number A is pre¬ 
dicted. However, often the calculations use simple forms 
for the matter distribution of the nucleus. 

In this Report we investigate the A-dependence us¬ 
ing realistic matter distributions. We will do this for 
two schematic models, covering the extremes sketched 
above. The hrst one (called model I) assumes a purely 
partonic mechanism, hadronization occurring outside of 
the nucleus (point B in Fig. ^ effectively at infinity) and 
thus absorption of the produced hadrons playing no role. 
The nuclear attenuation is then caused by the energy 
loss of the produced quark due to multiple scattering and 
gluon bremsstrahlung, which gives rise to a change of the 
hadron fragmentation function and thus to nuclear at¬ 
tenuation. Because of the Landau-Pomeranchuk-Migdal 
interference effect ^ this energy loss depends quadrat- 
ically on the length of matter traversed by the hit quark 
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(for the present discussion we neglect contributions from 
other quarks; see Refs. HQ for more details). In terms 
of the picture of Fig. n]it depends on the square of the 
density-averaged distance the parton travels within the 
nucleus from the point where it is created (point A in 
Fig. P). 

On the other hand our second schematic model (called 
model II) assumes that a possible attenuation is com¬ 
pletely due to absorption of the produced hadron, so 
nothing happens between the time the parton is produced 
and the hadron is formed (points A and B in Fig.^). In a 
Glauber approach this attenuation depends on the cross 
section for absorption of the hadron^ and the density- 
averaged distance the hadron travels through the nucleus 
after it has been formed (see Fig.^). Taking the absorp¬ 
tion cross section to be constant (and small enough that a 
linear approximation is sufficient), the latter determines 
the nuclear attenuation. 

This means we have to calculate the following two in¬ 
tegrals: 
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(model I) (1) 
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Here the exponential models the distribution of the for¬ 
mation distances If, and the matter densities pA are nor¬ 
malized to A. (By entering these quantities with cor¬ 
responding dynamical factors into the appropriate for¬ 
mula’s of the models, hadron production cross sections 
and from these values for the attenuation Ra can be 
calculated. However, here we are interested in the A- 
dependence (and moreover the used models are extreme 
and schematic). Under the assumption that the cross sec¬ 
tion is linear in (t'j) and (t//), which is a good first-order 
approximation, the A-dependence of these quantities car¬ 
ries over into the one of the cross sections.) 

It can easily be shown that for a nucleus with a mass 
density distribution described by one scale parameter, as 
in the case of a rectangular (as in the liquid-drop model) 
or a Gaussian distribution, the value of (tj) is propor¬ 
tional to the equivalent radius or rms radius squared, 
which in those cases leads to an dependence. In 
case of model H one finds for (t//) for a rectangular dis¬ 
tribution an A^/^ dependence when Lf = 0, and a larger 



FIG. 2: Values of (tf) as function of A for 2-parameter Fermi 
(2pF) matter distributions (solid curve) and for actual nuclei 
(circles). The dashed line is a power-law (A“) fit to the lat¬ 
ter, excluding the point for the deuteron. The dotted line in 
addition excludes the point for “^He in the fit. 


exponent'^ (e.g., about 0.55 for Lf = 4 fm) at larger Lf, 
and similar for a Gaussian distribution. 

However, neither a rectangular nor a Gaussian distri¬ 
bution is a good representation of the mass distribution 
of a real nucleus. Therefore we have evaluated {tj) and 
{tij) for a more realistic distribution, described by a 2- 
parameter Fermi (Saxon-Woods) form 

p^(r)=po/[l + e-(’'-^)/“] (3) 

with parameters po = 0.170 nucleons/fm^, a = 0.5 fm, 
and the value of c so as to give a nucleus with A nucleons. 
(This form gives a reasonably good global description of 
the mass distribution down to low values of A). The 
results are given in Tableland are shown in Figs. Inland 

m 

It can be seen from Fig. [3 that for a 2pF matter dis¬ 
tribution the A-dependence of {t'j) is rather steep. If one 
tries to describe it with the power law A“, it would re¬ 
quire a value of a well in excess of 2/3. Fig. 0 shows 
for the case of model H that the slope of the curve, in 
essence the value of a, increases when Lf increases. Val¬ 
ues range from about a = 0.40 for L/ = 0 fm to a = 0.60 
for Lf = A fm. 

Given these findings, and since it is known that the 
parameters for actual nuclei are slightly irregular due to, 
e.g., shell closures, it is interesting to see the behavior of 


® In our simple model we do not take into account coupled-channels 
processes. Possible effects of these are discussed in Ref. |^. 


^ This effect of taking into account a distribution of formation 
distances has been noted already in Ref. |^. 
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TABLE I: Values of (if) and {tii) for nuclei of mass number 
A, calculated with a 2pF distribution (see Eq. |2J with half 
radius c (second column) and po = 0.170 fm“® and a = 0.5 
fm fixed. The values of (tu) were calculated for five different 
values of Lf (columns 4-8). 



A 


FIG. 3: Values of {tu) as function of A for Lf — 0,2,4 fm 
(solid curves, top to bottom) and Lf = 1, 3 fm (dashed curves, 
top to bottom) for 2pF matter distributions and for actual 
nuclei (symbols) for Lf — 1,3 fm only. The dot-dashed lines 
represent (arbitrarily normalized) and behaviors, 

as indicated. 


(tj) and (t//) for real nuclei, since those are used in exper¬ 
iments. For that purpose we have used parameterizations 
[11113 of measured charge distributions (since the neu¬ 
tron distribution is very similar the error introduced by 
using the charge distribution instead of the matter dis¬ 
tribution is small, and irrelevant for the conclusions of 
the present study) for the nuclei ^H, '*He, ^®Si, 

4°Ca, 4®Ca, s-^Kr, “^Sn, i32xe and 2°®Pb. The results 
are also given in Table HD and shown as the symbols in 
Figs. 121 and El 


TABLE II: Calculated values of (tj) and (tj/) for various 
nuclei. The rms radii of the employed mass distributions are 
given in the second column. The values of (t/r) were calcu¬ 
lated for five different values of Lf (columns 4-8). 


nucleus 

^ rms 

[fm] 

(i?) 

ffm-"l 

{tu) 

[fm-2] 






Lf [fm] 






0.0 

1.0 

2.0 

3.0 

4.0 


2.13 

0.011 

0.068 

0.041 

0.027 

0.021 

0.017 

"He 

1.72 

0.050 

0.159 

0.095 

0.063 

0.047 

0.038 

12c 

2.45 

0.086 

0.218 

0.153 

0.111 

0.086 

0.070 

16 o 

2.73 

0.101 

0.238 

0.172 

0.127 

0.100 

0.082 


3.08 

0.173 

0.316 

0.239 

0.181 

0.145 

0.120 

"°Ca 

3.48 

0.216 

0.356 

0.277 

0.215 

0.174 

0.146 

«Ca 

3.47 

0.309 

0.428 

0.333 

0.259 

0.210 

0.175 

«"‘Kr 

4.25 

0.395 

0.493 

0.403 

0.325 

0.270 

0.230 

“®Sn 

4.67 

0.523 

0.571 

0.475 

0.390 

0.328 

0.281 

«2Xe 

4.83 

0.570 

0.598 

0.500 

0.413 

0.348 

0.300 

208pb 

5.51 

0.811 

0.719 

0.615 

0.519 

0.445 

0.387 



FIG. 4: Values of {tu) as function of A for Lf =0,2,4 fm 
for actual nuclei (symbols). The lines are results of power-law 
(A“) fits, excluding the point for the deuteron. 

A first observation is that while the results for A > 4 
scatter around the global curve, the results for are 
way down. This is due to being a rather dilute system. 
Furthermore, it is seen that the results for ^He lie slightly 
above the global curve, which is related to "‘He being a 
relatively dense system. Fitting the A-dependence of the 
results for the real nuclei including "‘He (but leaving out 
^H), one finds an exponent of 0.74 for {t]) (see Fig. |21) 
and values between 0.40 and 0.60 for {tu) depending on 
the value of Lj (see fig. El). When "‘He is not included in 
the fit, the values of a in case of {tu) change by less than 
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0.02, whereas in case of [t]) the value of a increases from 
0.74 to 0.79. Thus, in trying to extract an >l-dependence 
from experimental data it matters if one uses ^He as low¬ 
est A nucleus, or, e.g., 

Given these results it would in principle be possible 
to discriminate on account of the ^-dependence between 
the two extreme mechanisms that we have used here. 
However, in practice the process of hadronization in a nu¬ 
cleus most probably will be a combination of these mech¬ 
anisms, with possibly even different dependences on path 
lengths in the nucleus than employed here. But whatever 
the mechanism, the H-dependence will be an important 
ingredient. Our results show that in model calculations 
of the H-dependence for comparison with experimental 
data, global mass distributions, such as a rectangular 
distribution employed in the liquid-drop model, are not 


adequate. Instead, the use of experimentally established 
density functions for the nuclei actually used in the ex¬ 
periment is essential, where it even matters whether ^He 
with its relatively large density is included or not. 
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